DECLARE FUNCTION Ei! (X)
'
'This is a 1-layer drawdown and buildup program with afterflow that utilizes
'US oilfield units with time in minutes and uses the Ei-function as 2 PD.
'
CLS
PRINT
'PRINT "Ready to enter print control parameter IPRN and input parameters"
'PRINT "   IPRINT = 1 sends drwdwn and bld data to printer and bld to file"
'PRINT "   IPRINT = 0 sends only bld data to file C:\PRESS.PRN"
PRINT
'
'INPUT "IPRINT= ", IPRN

'
'C2= 4.395E-06 for time in minutes; AK=0.0227528 minutes for Oil/gas Type Curves (oilfield units)

'

'INPUT "Enter: T/F=", R
'INPUT "Enter: DELTA-minutes = [0.1 * 10,000/(T/F)] ", D           ' = 0.1 * 10,000/(T/F)
'INPUT "Enter: TIMOFF-minutes = [2.5 * 2.2*10^6/(T/F)] ", TOF#      ' enter this value 2.5 * 2.2*10^6/(T/F)
'INPUT "Enter: R^2/C2*ETA1 = [0.0227528] ", AK
'INPUT "BUILDUP TIME-minutes = [10,000] ", TBD
'INPUT "SKIN = [0] ", S  ' enter 0 (zero)
INPUT "Enter: T/F=", R
PRINT ""
PRINT " DELTA-minutes        = [0.1 * 10,000/(T/F)] "  ', D           ' = 0.1 * 10,000/(T/F)
PRINT " TIMOFF-minutes       = [2.5 * 2.2*10^6/(T/F)] "      ', TOF#      ' enter this value 2.5 * 2.2*10^6/(T/F)
PRINT " R^2/C2*ETA1          = [0.0227528] "'  , AK
PRINT " BUILDUP TIME-minutes = [10,000] " '     , TBD
PRINT " SKIN                 = [0] ", '   S  ' enter 0 (zero)
PRINT ""

' ----------- calculate input values -----------
D = .1 * 10000 / R
TOF# = 2.5 * 2200000! / R
AK = .0227528
TBD = 10000
S = 0
' -----------------------------------------------



' ------ set default values -------
IF D = 0 THEN
   D = .1 * 10000 / (R)             ' -- delta minutes
END IF

IF TOF# = 0 THEN
   TOF# = 2.5 * 2.2 * 10 ^ 6 / (R)    ' -- Timeoff
END IF

IF AK = 0 THEN
   AK = .0227528#                    ' -- Constant
END IF

IF TBD = 0 THEN
   TBD = 10000                       ' -- Build-up time
END IF

' ---------------------------------


'
'IF IPRN = 0 THEN GOTO START
'LPRINT "THIS IS A ONE-LAYER AFTERFLOW DRAWDOWN AND BUILDUP PROGRAM"
'LPRINT
'LPRINT
'LPRINT "R^2/C2*ETA1=  ", AK
'LPRINT
'LPRINT "DELT=  "; D
'LPRINT
'LPRINT "T/F=  "; R
'LPRINT
'LPRINT "SKIN=  "; S'
'LPRINT
'LPRINT "TIMOFF=  "; TOF#
'LPRINT
'LPRINT "BUILDUP TIME =", TBD
'LPRINT
'LPRINT
'
'Setup the format for calculational results in columnar form on printer
'
'LPRINT "  TIME,MIN          Q(t)/Q         DELP*F/Q      SHUT-IN, MIN      RECOVERY"
'LPRINT "____________     ____________    ____________    ____________    ____________"
'LPRINT
'
'Three dimensional constants, C1, C2, and C3, appear in this program in
'order to accommodate any system of units.  The constant C2 is part of the
'calculated value of the input AK for the argument of the Ei-function and
'is chosen for time units in minutes.  The constant C1 multiplies the Ei
'and equals 1/4*PI in consistent units whereas C3 multiplies the storage
'factor and equals 1 in consistent units.  US oilfield values appear below.
'
START:
'
'Creation of ASCII file "A:\AFBLD.PRN" to save results of calculations
'for pressure recovery above flowing pressure after shut in.
'
'OPEN "C:\PRESS.PRN" FOR OUTPUT AS #1
PRINT
PRINT "Ready to set up data file (8-characters) "
PRINT
INPUT "FILE TITLE: ", AS$
AFile$ = AS$ + ".TXT"

OPEN AFile$ FOR OUTPUT AS #1

PRINT #1, "TITLE:  ", AS$
PRINT #1, "T/F............. ", R
PRINT #1, "Delta-T......... ", D
PRINT #1, "Time Off........ ", TOF#
PRINT #1, "Build-up Time... ", TBD
PRINT #1, "SI MINUTES", " DELP*F/Q"
PRINT #1, " ": PRINT #1, " "
'
'Calculations start at this point
'
DIM T(2000), DT(2000), Q(2000)
C1 = 70.6
C3 = 256.456
'
'The next calculation converts the input value for T/F to C1*C3*(F/T)
'and divides the value of AK by 4 for the Ei arguments and sets an index
'value for setting up the buildup parameters at shut-in time
'
C4 = C1 * C3 / R
AK = AK / 4
INDEX = 1
'
'Calculations for q/qo and Delpress start here
'
Q(1) = 0
T(1) = D
DT(1) = D
I = 1
'
'Recycle on I label is "RECYCLE"
'
RECYCLE:
'
SUM1 = 0
SUM2 = DT(I)
SUM3 = 0
X = AK / DT(I)
UI = Ei(X) + 2 * S
CII = DT(I) + C4 * UI
K = 1
DO UNTIL K = I
'
X = AK / (T(I) - T(K - 1))
U = Ei(X)
X = AK / (T(I) - T(K))
DU = U - Ei(X)
CIK = DT(K) + C4 * DU
SUM1 = SUM1 + CIK * Q(K)
SUM2 = SUM2 + DT(K)
SUM3 = SUM3 + Q(K) * DT(K)
K = K + 1
LOOP
'
IF T(I) > TOF# THEN SUM2 = TOF#
Q(I) = (SUM2 - SUM1) / CII
SUM3 = SUM3 + Q(I) * DT(I)
DPI = (SUM2 - SUM3) / C3

IF T(I) > TOF# THEN GOTO BUILDPRN

'IF IPRN = 0 THEN GOTO FLOW
'LPRINT USING "##.####^^^^"; T(I);
'LPRINT "      "; USING "##.####^^^^"; Q(I);
'LPRINT "     "; USING "##.####^^^^"; DPI
'
FLOW:
'
Q = Q(I)
IF T(I) > TOF# THEN Q = 1 - Q
I = I + 1
IF I > 2000 THEN PRINT "Calculations exceeded the permitted 2000 time steps"
IF I > 2000 THEN GOTO FEND
DT(I) = D
IF Q > .25 THEN DT(I) = 5 * D
IF Q > .5 THEN DT(I) = 10 * D
IF Q > .75 THEN DT(I) = 20 * D
IF Q > .9 THEN DT(I) = 50 * D
IF Q > .96 THEN DT(I) = 100 * D
IF Q > .99 THEN DT(I) = 250 * D
T(I) = T(I - 1) + DT(I)
IF T(I) > TOF# + TBD THEN PRINT #1, " "
IF T(I) > TOF# + TBD THEN PRINT #1, "FLOW TIME =", TOF#
IF T(I) > TOF# + TBD THEN PRINT #1, "DELP*F/Q at shut-in=", DPIO
IF T(I) > TOF# + TBD THEN GOTO FEND
IF T(I) > TOF# THEN GOTO BUILDUP
GOTO RECYCLE
'
BUILDPRN:

'IF IPRN = 0 THEN GOTO SKIP
'LPRINT USING "##.####^^^^"; T(I);
'LPRINT "      "; USING "##.####^^^^"; Q(I);
'LPRINT "     "; USING "##.####^^^^"; DPI;
'LPRINT "     "; USING "##.####^^^^"; T(I) - TOF#;
'LPRINT "     "; USING "##.####^^^^"; DPIO - DPI
'
SKIP:
Data.Point = Data.Point + 1
LOCATE 18, 5
PRINT "Data Point # "; Data.Point
PRINT #1, T(I) - TOF#, DPIO - DPI
GOTO FLOW
'
BUILDUP:
'
'Setup for pressure buildup calculations done only one time at shut in
'
IF TBD = 0 THEN GOTO FEND
IF INDEX = 0 THEN GOTO RECYCLE
TOF# = T(I - 1)
DPIO = DPI
DT(I) = D
INDEX = 0
T(I) = T(I - 1) + DT(I)
GOTO RECYCLE
'
FEND:
PRINT #1, "Data Points "; Data.Point
LOCATE 22, 5
PRINT "Progam End <Press any key>"
DO
LOOP UNTIL INKEY$ <> ""

RESET
SYSTEM


END

FUNCTION Ei! (Z)
IF Z = 0 OR Z < 0 THEN E1 = 1E+38
IF Z = 0 OR Z < 0 THEN PRINT "EI UNDEFINED FOR ZERO OR NEGATIVE ARGUMENTS"
IF Z = 0 OR Z < 0 THEN GOTO FINI
IF Z = 1 OR Z > 1 THEN GOTO BIGZ
'
'USES SERIES APPROXIMATION WITH ERROR 2E-08
CONST C2 = .577215767#
I1 = 1: FAC = 1: T1 = 1: E1 = 0
DO UNTIL ABS(T1) < 2E-08
FAC = I1 * FAC
T1 = ((-Z) ^ I1) / (I1 * FAC)
E1 = E1 - T1
I1 = I1 + 1
LOOP
E1 = E1 - LOG(Z) - C2
GOTO FINI
'
BIGZ:
'USES HASTINGS, P.190,WITH ERROR 2E-08
CONST A0 = .2677737343#: CONST B0 = 3.9584969228#
CONST A1 = 8.634760892499999#: CONST B1 = 21.0996530827#
CONST A2 = 18.059016973#: CONST B2 = 25.6329561486#
CONST A3 = 8.573328740099999#: CONST B3 = 9.573322345399999#
'
E1 = A0 + A1 * Z + A2 * (Z ^ 2) + A3 * (Z ^ 3) + (Z ^ 4)
E1 = E1 / (B0 + B1 * Z + B2 * (Z ^ 2) + B3 * (Z ^ 3) + (Z ^ 4))
E1 = E1 * EXP(-Z) / Z
GOTO FINI
'
FINI:
'FUNCTION DEFINITION
Ei = E1
END FUNCTION

